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ABSTRACT 

We study minimal scenarios of resonant leptogenesis near the electroweak phase transition. These 

models offer a number of testable phenomenological signatures for low-energy experiments and 

Q_i' future high-energy colliders. Our study extends previous analyses of the relevant network of 

Q^' Boltzmann equations, consistently taking into account effects from out of equilibrium sphalerons 

\ and single lepton flavours. We show that the effects from single lepton flavours become very 

important in variants of resonant leptogenesis, where the observed baryon asymmetry in the 
• i-H , 

^ ■ Universe is created by lepton-to-baryon conversion of an individual lepton number, for example 
' 

\ that of the r-lepton. The predictions of such resonant r-leptogenesis models for the final baryon 
asymmetry are almost independent of the initial lepton-number and heavy neutrino abundances. 
These models accommodate the current neutrino data and have a number of testable phenomeno- 
logical implications. They contain electroweak-scale heavy Majorana neutrinos with appreciable 
couplings to electrons and muons, which can be probed at future e~^e~ and high-energy 
colliders. In particular, resonant r-leptogenesis models predict sizeable Ou(3(3 decay, as well as e- 
and yU-number-violating processes, such as /x — > 67 and ;U — > e conversion in nuclei, with rates 
that are within reach of the experiments proposed by the MEG and MECO collaborations. 

PACS numbers: 11.30.Er, 14.60.St, 98.80.Cq 
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1 Introduction 



The origin of the baryon asymmetry in our Universe (BAU) has always been one of 
the central topics in particle cosmology. Recently, the high-precision determination of 
many cosmological parameters, including the baryon-to-photon ratio of number densities, 
?7b ~ 6.1 X 10^^° [1], has given renewed momentum for extensive studies on this topic [2]. 
The established BAU provides one of the strongest pieces of evidence towards physics be- 
yond the Standard Model (SM). One interesting suggestion for explaining the BAU, known 
as leptogenesis [3], is linked with neutrinos. Although strictly massless in the SM, neutrinos 
can naturally acquire their small observed mass through the presence of superheavy part- 
ners and the so-called seesaw mechanism [4] . These superheavy neutrinos are singlets under 
the SM gauge group and may therefore possess large Majorana masses that violate lepton 
number (L) conservation by two units. In an expanding Universe, these heavy Majorana 
neutrinos will in general decay out of equilibrium, potentially generating a net lepton asym- 
metry. The so-produced lepton asymmetry will eventually be converted into the observed 
BAU [3] by means of in-thermal equilibrium [B + L)-violating sphaleron interactions [5]. 

One difficulty faced by ordinary seesaw models embedded in grand unified theories 
(GUTs) is associated with the natural mass scale of the heavy Majorana neutrinos. This is 
expected to be of order the GUT scale Mqut = 10^^ GeV. On the other hand, inflationary 
supergravity models generically predict a reheating temperature T^ch oi order 10'^ GeV. 
In these models, a significant constraint on the upper bound for Tj-eh comes from the 
requirement that gravitinos are underabundant in the early Universe and so their late decays 
do not disrupt the nucleosynthesis of the light elements [6]. However, the low Treh gives 
rise to another constraint within the context of thermal leptogenesis. The heavy Majorana 
neutrino, whose L- violating decays are responsible for the BAU, has to be somewhat lighter 
than Treh ~ 10^ GeV, so as to be abundantly produced in the early Universe. Such a mass 
for the heavy Majorana neutrino should be regarded as unnaturally low for GUT-scale 
thermal leptogenesis. Finally, further constraints on successful GUT-scale leptogenesis [7- 
10] may be obtained from solar and atmospheric neutrino data [11]. 

The aforementioned problem with a low reheating temperature may be completely 
avoided in models that reahze low-scale thermal leptogenesis [12-14]. In particular, the 
lowering of the scale may rely on a dynamical mechanism, in which heavy-neutrino self- 
energy effects [15] on the leptonic asymmetry become dominant [16] and get resonantly 
enhanced [12], when a pair of heavy Majorana neutrinos has a mass difference comparable 
to the heavy neutrino decay widths. In [17], this dynamical mechanism was termed resonant 
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leptogenesis (RL). As a consequence of RL, the heavy Majorana mass scale can be as low as 
~ 1 TeV [12, 13] in complete agreement with the solar and atmospheric neutrino data [17]. 

A crucial model-building aspect of RL models is that such models have to rely on a 
nearly degenerate heavy neutrino mass spectrum. Although, without any additional lepton- 
flavour symmetry, such a requirement would appear very fine-tuned, there is no theoretical 
or phenomenologically compelling reason that would prevent the singlet neutrino sector of 
the SM from possessing such a symmetry. Specifically, the RL model discussed in [12], which 
was motivated by Eg unified theories [18], was based on a particular lepton symmetry in 
the heavy neutrino sector. This lepton symmetry was broken very approximately by GUT- 
and/or Planck-scale suppressed operators of dimension 5 and higher. In [17], another RL 
scenario was put forward based on the Froggatt-Nielsen (FN) mechanism [19], where two 
of the heavy neutrinos naturally had a mass difference comparable to their decay widths. 
Recently, several constructions of RL models appeared in the literature within the context 
of supersymmetric theories [20-23], or even embedded in SO(IO) unified theories [24,25]. 

One of the great advantages of RL models is that their predictions for the BAU 
are almost independent of the primordial L-number, S-number and heavy neutrino abun- 
dances [13,17]. This fact may be explained as follows: in RL scenarios, the L- violating 
decay widths of the heavy Majorana neutrinos can be significantly larger than the Hubble 
expansion rate H of the Universe. As a consequence, the heavy Majorana neutrinos can 
rapidly thermalize from their decays, inverse decays and scatterings with the other SM par- 
ticles in the plasma, even if there were no heavy Majorana neutrinos at high temperatures. 
Moreover, in this high temperature regime, any pre-existing lepton asymmetry will rapidly 
be driven to zero, due to the L- violating inverse decays and scattering processes which 
are almost in thermal equilibrium. As the Universe cools down, a net lepton asymmetry 
can be created at temperatures just below the heavy neutrino mass as a consequence of 
the aforementioned CP-violating resonant enhancement that occurs in RL models. This 
L asymmetry will survive wash-out effects and will be converted by the {B + L)-violating 
sphalerons into the observed BAU. 

In this paper we provide a detailed study of a new variant of RL where a given 
single lepton fiavour asymmetry is resonantly produced by the quasi-in-equilibrium decays 
of heavy Majorana neutrinos of a particular family type. Such a variant of RL was first 
discussed in [26], and for the case of the r-lepton number this mechanism has been called 
resonant t -leptogenesis (RrL). This mechanism makes use of the property that, in addition 
to B — L, sphalerons preserve the individual quantum numbers ^B—Lg^^T- [27-30]. In a RrL 
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model, the generated excess in the Lr number will be converted into the observed BAU, 
provided the L,--violating reactions are not strong enough to wash out such an excess. 

Although our focus will be on minimal non- super symmetric 3-generation RL mod- 
els, supersymmetry could account for the origin of the electroweak-scale heavy Majorana 
neutrinos. In particular, one may tie the singlet Majorana neutrino mass scale mjsf to the 
//-parameter through the vacuum expectation value (VEV) of a chiral singlet superfield 
S [31]. The proposed model is a variant of the so-called Next-to- Minimal Supersymmetric 
Standard Model (NMSSM) and is described by the following superpotential (summation 
over repeated indices implied): 

W = Wmssm{i^^^) + K^U%VjR + XS%H2 + ^SViRViR + (1.1) 

where VrMSSM(Ai = 0) is the superpotential of the well-known Minimal Supersymmetric 
Standard Model (MSSM) without the //-term, and ^1,2,3 and z?i,2,3r are the Higgs- 
doublet, lepton-doublet and right-handed neutrino superfields, respectively. Once the scalar 
component of S develops a VEV vs, then both the would-be //-parameter, /i = Af 5, and the 
S0(3)-symmetric singlet scale, rriN = ^pvs, are expected to be comparable in magnitude 
(asumming that A ~ p), thus providing a natural framework for the possible existence 
of 3 nearly degenerate electroweak-scale heavy Majorana neutrinos [32]. In this minimal 
extension of the MSSM, the predictions for the BAU will depend on the size of the soft 
SUSY-breaking mass scale Msusy- However, if Msusy is relatively larger than the singlet 
Majorana neutrino mass scale mjv, e.g. Msusy ~ 2miv, the dominant source of leptogenesis 
will be the minimal non-supersymmetric sector that we are studying here, so our predictions 
will remain almost unaffected in this case. 

As mentioned above, single lepton-flavour effects on the net L and B asymmetries play 
a key role in RrL models. To properly treat these as well as SM chemical potential effects, 
the relevant network of the Boltzmann equations (BEs) needs to be extended consistently. 
In particular, single lepton-flavour effects can have a dramatic impact on the predictions for 
the B asymmetry. These predictions for the BAU can differ by many orders of magnitude 
with respect to those obtained in the conventional BE formalism, which is commonly used in 
the literature. Although our primary interest will be to analyze RL models, we should stress 
that single lepton-flavour effects could also signiflcantly affect the predictions obtained in 
hierarchical leptogenesis scenarios. The improved set of BEs derived here will therefore be 
of general use. 

Another important question we wish to address is whether the leptogenesis scale 
can be lowered to energies 100-250 GeV, very close to the critical temperature Tc, where 
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the electroweak phase transition occurs. In this temperature region, freeze-out effects 
from sphaleron processes dropping out of equihbrium need to be considered, as they can 
significantly modify the predicted values for the final baryon asymmetry. Our treatment of 
these sphaleron freeze-out effects will be approximate and based on the calculations of [27, 
30,33]. Our approximate treatment is motivated by the fact that, within the framework 
of RL models, the creation of a net lepton asymmetry at the electroweak scale does not 
require the electroweak phase transition to be strongly first order. 

Most importantly, in models where the BAU is generated from an individual lepton- 
number asymmetry, a range of testable phenomenological implications may arise. The key 
aspect is that scenarios such as RrL can contain heavy Majorana neutrinos with appreciable 
Yukawa couplings to electrons and muons. The (normalized to the SM) ly^-boson couplings 
of e and /x leptons to these heavy Majorana neutrinos could be as large as 10^^. For 
electroweak- scale heavy neutrinos, such couplings would be sufficient to produce these 
particles at future e^e~ and iJi^ ijT colliders. Furthermore, minimal (non-supersymmetric) 
3-generation RrL models can predict /i — >■ 67 and fi ^ e conversion in nuclei at rates that 
can be tested by the foreseeable experiments MEG at PSI [34] and MECO at BNL [35], 
respectively. Finally, RrL models naturally realize an inverted hierarchy for the light 
neutrino spectrum and therefore also predict neutrinoless double beta {Oi>f3f3) decay with 
a sizeable effective neutrino mass as large as ~ 0.02 eV. This value falls within reach 

of proposals for future Oz//?/5-decay experiments sensitive to \{'m)\ ~ 0.01-0.05 eV [36], 
e.g. CUORE (i30Te), GERDA (^^Ge), EXO (^^exe), MOON (^o^Xe), XMASS (^^^Xe), 
CANDLES (^^Ca), SuperNEMO (^^Se) etc. 

Our paper has been organized as follows: Section 2 presents a minimal model for 
resonant r-leptogenesis. In Section 3 we derive the EEs for single lepton flavours, by 
carefully taking into account SM chemical potential effects. Technical details pertinent to 
this derivation have been relegated to Appendix A. In Section 4 we review the calculation 
of out of equilibrium sphaleron effects at the electroweak phase transition and apply it 
to leptogenesis. In Section 5 we give several numerical examples of RrL models, focusing 
our attention on scenarios that can be tested at future e'^e~ and colliders and in 

low-energy experiments. In particular, in Section 6, we present predictions for lepton- 
flavour-violating (LEV) processes, such as — > 67, /i — > eee and — > e conversion in 
nuclei. Finally, we present our conclusions and future prospects in Section 7. 
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2 Minimal Model for Resonant r-Leptogenesis 



There have been several studies on RL models in the literature [13,14,17,20,21,23,25]. 
Here, we will focus our attention on a variant of resonant Icptogenesis where the BAU 
is generated by the production of an individual lepton number [26]. For definiteness, we 
consider a minimal (non-supersymmetric) model for RrL. 

Let us start our discussion by briefly reviewing the relevant low-energy structure of 
the SM symmetrically extended with one singlet neutrino Uir per i family (with i = 1,2,3). 
The leptonic Yukawa and Majorana sectors of such a model are given by the Lagrangian 



1 ^ 

-^M,Y = o E ( i^iR)^ iMs)ij UjR + h.C. 



3=1 



3 



(2.1) 



i=e,n,T 

where Li = {un, liLy are the left-handed lepton doublets *, liR are the right-handed leptons, 
and $ is the isospin conjugate of the Higgs doublet 

In the Lagrangian (2.1), we have defined the individual lepton numbers ive,^,T in 
the would-be charged-lepton mass basis, where the charged-lepton Yukawa matrix is 
positive and diagonal. In fact, without loss of generality, it can be shown that sphaleron 
transitions exhibit a U(3) flavour symmetry and so they can be rotated to become flavour 
diagonal in the same would-be mass basis. To prove this, one may write the operator Ob^l 
responsible for 5 + L- violating sphaleron transitions as follows (group-invariant contraction 
of the colour and weak degrees of freedom implied) [28] : 

Ob+l = i{Q[Q[Q[L',, (2.2) 

i=l 

where Q[ and L'^ denote the quark and lepton doublets defined in an arbitrary weak basis. 
The operator Ob+l contains the fully antisymmetric operator combinations: Q'lQ^Q'z 3^nd 
L'^^L^L'^, which are invariant under U(3) fiavour rotations [37]. Thus, we can use this U(3)- 
rotational freedom to render the charged lepton and up-quark sectors fiavour diagonal and 
positive. 

To obtain a phcnomcnologically relevant model, at least 3 singlet heavy Majorana 
neutrinos 1^1^2,3 r ai'c needed and these have to be nearly degenerate in mass. To ensure the 
* Occasionally wc will also denote the individual lepton numbers with Lp ^ t-. This apparent abuse of 
notation should cause no confusion to the reader, as the precise meaning of L^^^^r can be easily inferred 
from the context. 
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latter, we assume that to leading order, the heavy neutrino sector is SO (3) symmetric, i.e. 



Ms = thnU + ^Ms , (2.3) 

where I3 is the 3x3 identity matrix and AMs is a general S0(3)-breaking matrix. As 
we will discuss below, compatibility with the observed light neutrino masses and mixings 
requires that {AMs)ij/mi\f < 10^''. for electroweak-scale heavy Majorana neutrinos, i.e. for 
m^v ~ 0.1-1 TeV. One could imagine that the soft S0(3)-breaking matrix AMs originates 
from a sort of Proggatt-Nielsen mechanism [19]. 

In order to account for the smallness of the light neutrino masses, we require that the 
neutrino Yukawa sector possesses a leptonic symmetry. This will explicitly break the 
imposed S0(3) symmetry of the heavy neutrino sector to a particular subgroup S0(2) ~ 
U(l)i. For example, one possibihty relevant to RrL is to couple all lepton doublets to a 
particular heavy neutrino combination: {u2R + i^zR)- In detail, the U(l)i charges of the 
fields are 

Q{U) - Q{hn)-1, - -^(^^^^^) - 1' Qi^.n)-0. (2.4) 

As a result of the U(l)/ symmetry, the matrix for the neutrino Yukawa couplings reads: 

^0 ce-*'^/'' ce'''/'^ J 



h 



+ . (2.5) 



In the above, a, h and c arc arbitrary complex parameters of the model. For electroweak- 
scale heavy neutrinos, the absolute value of these parameters has to be smaller than 
about 10~^, for phenomenological reasons to be discussed below and in Section 6. In 
particular, the requirement that an excess in Lr is protected from wash-out effects leads to 
the relatively stronger constraint |c| ^ 10~^. In addition, 5h^^ is a 3 x 3 matrix that pa- 
rameterizes possible violations of the U(l)z symmetry. It should be noted that the charged 
lepton sector and the leading S0(3)-invariant form of the heavy neutrino mass matrix are 
compatible with the U(l); symmetry. 

In this paper we shall not address the possible origin of the smallness of the S0(3)- 
and U(l) /-breaking parameters {AMs)ij and 5h\j-, as there are many different possibihties 
that could be considered for this purpose, e.g. the Froggatt-Nielsen mechanism [17,19], 
Planck- or GUT-scale lepton-number breaking [12,18]. Instead, in our model-building we 
will require that the symmetry breaking terms do not induce radiative effects much larger 
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than the tree-level contributions. This naturalness condition will be applied to the light 
and heavy neutrino mass matrices and Ms, respectively. 

We start by observing that the U(l); symmetry is sufficient to ensure the vanishing of 
the light neutrino mass matrix [38]. In fact, if U(l); is an exact symmetry of the theory, 
the light neutrino mass matrix will vanish to all orders in perturbation theory [39,40]. To 
leading order in the U(l)/-breaking parameters AM5, the tree-level light neutrino mass 
matrix m'^ is given by 



2 2 
2 2m N 



rriN 



(2.6) 



where v — 2Mw/gw = 245 GeV is the vacuum expectation value of the SM Higgs field 
As a minimal departure from U(l); in the neutrino Yukawa sector, we consider that this 
leptonic symmetry is broken only by through 



I 



\ 



o^ 


Oy 



(2.7) 



In this case, the tree-level light neutrino mass matrix (2.6) takes on the form 



m 



2mjv 



1jVq2 



V 



ab — £eS,, 

ruN 



TUN 

AiUN ip. _ g.2 

^bc- e.Sr 



ac - EeEr ^ 

ruN 

be - e^Er 

TUN ^ 



Am^v ^2 _ g,2 

TOJV 



(2.8) 



where Am^r = 2(AM5)23 + i[{AMs)33 — (AM5)22]- It is interesting to notice that in 
this type of 11(1); breaking, the parameters £e,^t,r enter the tree-level light neutrino mass 
matrix m" quadratically. As a consequence, one finds that for itln ~ v, these U(l)/- 
breaking parameters need not be much smaller than the electron Yukawa coupling he ~ 
10^^. Moreover, one should observe that only a particular combination of soft S0(3)- and 
U(l)rbreaking terms {AMs)ij appears in m'^ through Amjv. Nevertheless, for electroweak- 
scale heavy neutrinos with mass differences \AmN\/mi\f < 10^^, one should have \a\, \b\ < 
10^^ to avoid getting too large light neutrino masses much above 0.5 eV. As we will see 
more explicitly in Section 5, for the RrL scenario under study, the favoured solution will 
be an inverted hierarchical neutrino mass spectrum with large Ue^/j. and f^i/,- mixings [11]. 

In addition to the tree level contributions given in (2.8), there are Z- and Higgs-boson- 
mediated threshold contributions Sm" to m'' [39]. The contributing graphs are displayed 
in Fig. 1. In the heavy neutrino mass basis, where Ms = diag(miVi , ttin^ > "^ats), with 
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iVi,2,3 

(a) 



H 



/ \ 
I \ 

_1 L 



(b) 



Figure 1: Finite radiative effects contributing to the light-neutrino mass matrix. 



iTT'Ni ^ ^ ^Na, and h^'^ = h'^, these finite radiative corrections may conveniently be 
expressed as follows [40]: 



327r 



E 

a=l,2,3 



7)2 



Bo{0,m'M'z) - i?o(0,0,Mi) 



Ml 



(2.9) 



where — gl,/{ATT) and Mh is the SM Higgs boson mass. In (2.9), i?o(0,m^,m|) is the 
usual Pasarino-Veltman one-loop function [41], i.e. 



S„(0,m?,mS = Cw + 1 - Inf^l + In f !!!h) , 



\ fj? J m\ — m2 \miJ 



(2.10) 



and Cuv is a UV infinite constant. Moreover, in writing (2.9), we have neglected terms of 
order [(/i^q,)^ '^^]/"^lf„ > which are suppressed by higher powers of the small Yukawa couplings. 
It can easily be verified that the radiative lepton-number- violating contribution 5m^ to the 
light neutrino mass matrix is UV finite and /x-scale independent. For ^ and 
("^Ar„ — JTijvJ/mATi <^ 1, the expression (2.9) evaluated in the original weak basis simplifies 
to 

;2 h'^^AMsih'^^f 



3M| 



167r M'yY rriM 

For electroweak-scale heavy Majorana neutrinos m^^ ^ v and = 120-200 GcV, one 
may estimate that for (mjv„ — mN^)/m]^^ < lO^'' and |a|, \b\ < 10~^, the finite radiative 
effects 5m^ stay well below 0.01 eV. In fact, up to an overall coupling-suppressed constant, 
these corrections have the same analytic form as the first term on the RHS of (2.6). They 
can be absorbed by appropriately rescaling Am at defined after (2.8). As a consequence, 
these finite radiative effects do not modify the parametric dependence of the tree-level light 
neutrino mass matrix given in (2.8). 



(2.11) 
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We now turn our attention to the heavy Majorana neutrino sector. In this case, 
renormahzation-group (RG) running effects [23,42] become very significant. These effects 
exphcitly break the S0(3)-symmetric form of the heavy neutrino mass matrix, Ms{Mx) = 
mi\f I3, imposed at some high energy scale Mx, e.g. at the GUT scale. A fairly good 
quantitative estimate of these S0(3)-breaking effects can be obtained by solving the RG 
equation for the heavy neutrino mass matrix M5: 



dMs 
IT 



1 



167r2 



(2.12) 



with t — ln(Mx///). Considering that h^^ has only a mild RG-scale dependence and 
assuming that Ms{Mx) = mjv I3 at some high scale Mx, we may calculate the RG effects 
by running from Mx to the low-energy scale uin ^ v through the relation 



M5(mjv) - 

= Ms{Mx) - 



Ms{Mx) - ?^Re 



87r2 

|aP + |6|2 



87r2 



tun In 



Mx 
rriN 



In 



diag (0,1,1) + O 



(|a|2 + |6|2)i/2 



(2.13) 



If the scale Mx of the SO (3) symmetry imposed on Ms{Mx) is to be naturally associated 
with the scale Mqut ~ 10^^ GeV relevant to GUT dynamics, it can be estimated from (2.13) 
that the mass splittings \mi\f.^ — mpf^l/mj^]- and jmATg — mArJ/mjv should be larger than 10~^ 
for \a\, \b\ ~ 10~^ (|c|, |£e,/i,T| ~ 10~^). Instead, the mass difference \mp^^—mN^\/misf should 
be comparatively much smaller, as it is protected by an approximate U(l)i symmetry. In 
particular, we find that Itjin^ — TnN2\I^N — C'd^c/i.ral, |£e,/i ,rh\) < 10-'^. At this point we 
should stress that in the scenarios we are considering, RG effects predominantly modify 
the entries [AMs)u (with i = 1, 2, 3) in (2.3) and so they do not affect the fight neutrino 
mass matrix (2.8). However, these effects may affect the single lepton flavour asymmetries 
and the flavour-dependent wash-out factors that are discussed in the next section. 

In addition to RG effects, one might worry that thermal effects could significantly 
modify the heavy neutrino mass spectrum. However, thermal effects respect the underlying 
symmetries of the theory, such as global, chiral and gauge symmetries [43]. Hence, their 
impact on the heavy neutrino mass spectrum is controlled by the size of the SO (3)- and 
U(l)rbreaking parameters in the Yukawa neutrino sector. In the hard thermal loop (HTL) 
approximation [43,44], thermal corrections give rise to an effective heavy neutrino mass 
matrix Ms{T), which differs from the one evaluated at T = by an amount [17] 



Ms{T) - Ms(0) 



16 



Re 



rriN 



(2.14) 
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By comparing (2.14) with (2.13), we notice that thermal corrections have a parametric 
dependence very similar to the RG effects and are opposite in sign. Nevertheless, if the 
S0(3)-breaking scale Mx is identified with Mgut? RG effects become larger than thermal 
effects by at least a factor 3, for the temperature regime relevant to leptogenesis T < rriN. 

In Section 5 we will present numerical estimates of the BAU for electroweak-scale RL 
models that are motivated by the naturalness of the light and heavy neutrino sectors. As 
we mentioned above, this condition provides a potential link between these models and 
GUT-scale physics. 

3 Boltzmann Equations for Single Lepton Flavours 

In this section we derive a set of coupled BEs for the abundances of heavy Majorana 
neutrinos and each lepton fiavour. We follow a procedure analogous to the one presented 
in [17], where a number of controllable approximations were made. In particular, we assume 
Maxwell-Boltzmann statistics for the heavy Majorana neutrinos. For the SM particles, we 
instead consider the proper Bose-Einstein and Fermi-Dirac statistics, but ignore conden- 
sate effects [45]. The above simplifications are expected to introduce an error no larger 
than 20%. Furthermore, we neglect thermal effects on the collision terms, which become 
less significant in the temperature domain T < mjvj relevant to RL. As we will see more 
explicitly in Section 5, the latter approximation may partially be justified by the obser- 
vation [17] that the resulting BAU predicted in RL models is independent of the initial 
abundances of the heavy neutrinos and any initial baryon or lepton asymmetry. 

Various definitions and notations that will be useful in deriving the BEs are introduced 
in Appendix A. Adopting the formalism of [45,46], the evolution of the number density, 
Ua, of all particle species a can be modelled by a set of BEs. These are coupled first order 
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differential equations and may be generically written down as 



+ ZHug ^ - Y: 7(a^' ^Y) - ^^{Y^ aX') 



. . . . (3-1) 

where all possible reactions of the form aX' — > y or y — > aX' , in which a can be created 
or annihilated are summed over. If a is unstable, it could occur as a real intermediate state 
(RIS) in a resonant process like X — > a — > y. In this case, special treatment is required to 
avoid overcounting processes. 

In principle, there is a large number of coupled BEs, one for each particle degree 
of freedom. This number can be drastically reduced by noting that rapidly occurring 
SM processes hold most of the different particle degrees of freedom and particle species 
in thermal equilibrium. The non-zero chemical potentials of the particle species other 
than heavy Majorana neutrinos and leptons produce effects of 0{1) on the final baryon 
asymmetry [48]. These effects will be consistently included in the BEs for the heavy 
Majorana neutrinos A^i,2,3 and the lepton doublets i^e,/i,r- 

Although an infinite scries of collision terms could be added to each BE, only a few 
will have a significant contribution. Following the procedure in [17], terms of order hJ^-^h^ 
and higher will be neglected, where /ij. ~ h'^ arc the one loop resummcd effective Yukawa 
couplings introduced in [17]. Also neglected are terms of order It^ for 2 •e^ 2 scatterings 
with two external heavy Majorana neutrinos. This leaves 1^2 decays and inverse decays 
of heavy Majorana neutrinos and 2 2 scatterings between heavy Majorana 

neutrinos, lepton doublets, gauge bosons, quarks and the Higgs field, which are formally of 
order h^^g^, h^^g'^ and {h^fhl. 

An important step in the following derivation is the proper subtraction of RISs. For 

example, the process <-> L^^^ will contain real intermediate heavy Majorana neutrino 

^This formalism neglects coherent time-oscillatory terms describing particle oscillations in terms of num- 
ber densities, as well as off-diagonal number densities n^^, for the destruction of a particle species h and 
the correlated creation of a particle species a, where a and b could represent the 3 lepton flavours or the 
3 heavy neutrinos Ni^2,z- Although these effects can be modelled as well [47], their impact on the BAU 
is expected to be negligible. Specifically, coherent time-oscillatory terms between heavy Majorana neutri- 
nos will rapidly undergo strong damping, as a consequence of the quasi-in-thermal equilibrium dynamics 
governing RL models. This results from the fact that the decay widths TjVi^s.s of the heavy neutrinos are 
much larger than the expansion rate of the Universe. Additionally, the correlated off-diagonal number den- 
sities n^i will bo Yukawa- coupling suppressed 0{{h}')'^) with respect to the diagonal ones na,b, if the heavy 
neutrinos and the charged leptons arc defined in the diagonal mass basis. In particular, the contribution 
of n^f, to na,6 will be further suppressed 0{{h"Y). We will therefore neglect the effects of the coherent 
time-oscillatory terms and the off-diagonal number densities n^g on the BEs. 
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states. Their inverse decay and subsequent decay have aheady been accounted for in the 
BEs and must be subtracted to ensure that unitarity and CPT are respected [45]. 

In analogy to 2 <-> 2 scatterings, 2 — > 3 processes, such as LjQ'^ L^$^ii'-^, may 
also contain the heavy neutrinos Na as RISs. The resonant part of such a process consists 
of the reaction LjQ^ Nau'~"', followed by the decay — > L'j^^'^ . As before, to avoid 
double counting, we subtract the RISs from such a 2 ^ 3 process. Although the off-shell 
2-^3 process is a higher order effect than those we are considering, the subtracted resonant 
part contributes terms of order (h'^)'^hf^ and must be consistently included within the given 
approximations for the BEs. Specifically, the following relations among the collision terms 
are derived: 

a=l l = e,iJ,,T 

3 



a=l l = e,fj,,T 

a=l l = e,a,T 



2 

a=l l = e,ix,T 

(3.2) 

where a prime denotes subtraction of RISs, the indices j, k — e, ii,t label lepton flavour, 
and Sjk = (1 + Sjk)"^ is a statistical factor that corrects for the production or annihilation 
of identical lepton flavours. In addition, we have deflned the individual lepton-flavour 
asymmetries and branching ratios as 



Bat — 



r(Ar,^L;<l>)+r(iV,^Lf$t) 



'N^ - 7r^::77: t^-^^tt; ttttttt • (3-3) 



As CP violation in these processes is predominantly caused by the resonant exchange of 
heavy Majorana neutrinos, the CP-violating collision terms have been approximated in 
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terms of the CP-conserving ones as follows: 

l = e,fi,T l = e,iJ,,T 

K:"' = E^q:c' etc. (3.4) 

/ = e,/i,T 

Unlike the 2 — > 3 reactions, 3 — > 2 processes are treated differently. Although 3 — > 2 
processes could contain real intermediate A^q, states, coUision terms for their associated 
annihilation processes have not been included before. For example, in the process LjL^^ — > 
Qu'^, a real intermediate state could be coherently created from L and $ states. This 
coherent RIS would then interact with another L state producing Q and u'-^. Previously, 
the process N^L — > Qu'-^ has only been considered for heavy neutrinos in a thermally 
incoherent state. This implies that 3 — > 2 processes containing as RISs have not yet 
been accounted for and should not be subtracted. With the help of CPT and unitarity, 
one may therefore obtain the following relations for the 3^2 processes: 

^{L^^u ^ L^Q) - 7(i^f ^^^^^ - L^Q"") = Oih'^'hl) , 

7(L,$i.^ ^ L,g^) - 7(Lf ^ L^g) = 0{h^'hl). (3.5) 

As a consequence of this, 3 — > 2 processes will contribute extra CP-conserving 2—^2 
collision terms, through the resonant exchange of real intermediate A^^ states. Applying 
the narrow width approximation, we find 



a;=l 1 = e,fj,,T 

o=l l = e,ii,T 



c , 



e,tJ,,T 



a=l l = e,fj,,T 



c 

a=l l = e,ix,T 

(3.6) 

We may now employ (3.1) and write down the BEs in terms of the number densities 
of heavy Majorana neutrinos njv„ and the lepton-doublet asymmetries riALe.^.r) 



dt 



+ 3i/nAr„ = M - j E (7l;* + 7q„c'= + 7l:qc + 7l,„ 
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iV„y^ N^Lk JV„<J.t 



cq 



(3.7) 



driAL 



a=l 



n 



eq 



Hi;.- 



+ 57z 



«=i 



E ( T^'i + 27^- + 27£f + 27^^^'^ + 27^ 



+ ^7q„c + ^7$ty^ + ^eq ( 7q„c + 7$ty^ ) ) 



+ E 

fe = e,/i,T 



7Lp$t + 7$t$t + 7l,* + 7$$t 



k = e,ix,T Ik 

3 



7^4+ + 7$t$t - 7l,.| 



LkL'? 



+ 27":c^ 



+ Z^'^Nc^'^Nc, 2^ I 7l;QC + l^Liu + 7L($ + 7l, 
l = e,iJ,,T 



^7$tv^ 



(3.8) 



In the above set of BEs, we have only kept terms to leading order in riALj/n^^, and imple- 
mented the relations given in (3.2)-(3.6). 

All SM species in the thermal bath, including the lepton doublets Lg.^^^, possess 
non-zero chemical potentials. These chemical potentials can be expressed in terms of the 
lepton-doublet chemical potentials only, under the assumption that SM processes are in 
full thermal equilibrium [28]. This analysis yields the following relations: 



//y = 0, 



4 15 

I = e,iJ,,T I = e,iJ,,T I = e,/i,T 



21 



E f^Li , 



(3.9) 



--e,iJ,,T 



where /ix denotes the chemical potential of a particle species x. The relations (3.9) can 
be used to implement the effects of the SM chemical potentials in the BEs. They result 
in corrections to the so-called wash-out terms in both the lepton and heavy neutrino BEs. 
At this point we should also note that the BEs in their present form are most accurate 
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above Tc. As T approaches Tc, the assumption that the sphaleron processes are in thermal 
equihbrium becomes less valid. This will result in 0{v/T) corrections to the relations in 
(3.9). The inclusion of the bulk of these corrections will be considered in the next section. 

To numerically solve the BEs, it proves convenient to introduce a number of new 
variables. In the radiation dominated epoch of the Universe relevant to baryogenesis, the 
cosmic time t is related to the temperature T through 

* = 2H^) ' (^■"') 

where ^ 

^ = ^ , H{z) « 17.2 X , (3.11) 

T Mpianck 

with Mpianck = 1.2 X 10^^ GeV. Also, we normalize the number density of a particle species, 
^a, to the number density of photons, n-y, thereby defining the new parameter 77a, 

with 

2T3 2mi 1 
n^\z) = = ^ J . 3.13 

TT^ TT^ Z-^ 

To allow the BEs to be written in a slightly more compact form, we will use the summation 
conventions 

7fy^ = E 7^y^ , r^AL = E ^al, , (3.14) 

l = e,ix,T l = e,ix,T 

where X and Y stand for any particle state other than Li and Na- 

Using (3.7)-(3.14) and incorporating corrections due to the SM chemical potentials, 
the BEs for heavy Majorana neutrinos and lepton doublets are written down 

dVNo, Z I VNo. \ frD{ak) , pS(aA;) p5 (afc) 

Zeq~ I [ ^ Yukawa + Gauge 



dz H{z 
2 



s:k f^D{ak) , f\S{ak) f^S (ak) 

'/ALfc Ojv„ I J- -r J- Yukawa Gauge 

k = e,ix,T 



(3.15) 



di]. 



ALi 



dz 



H{z 
2 



-, \ ) 2^ "N^ \ cq -L I J- + i Yukawa + ^ 

J- J I a=l V^JVc / k = e,ti,T ^ 



S (ak) 
Gauge 



L a=l 



Er>3 I r^(aj) I p5(aj) (aj) (ai) , pl4^ (aj) 

-°iV„ \ ^ "I" ^ Yukawa ^ Gauge ^ Yukawa ^ Gauge 
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where 



{al) 

p£) {al) 

p5 (al) 
Yukawa 

pS (al) 
Yukawa 



(al) 
Yukawa 



-.S (al) 
Gauge 



^5 (al) 
Gauge 



^5 (al) 
Gauge 



-,W (at) 
Yukawa 



+ E 

k = e,fj,,T 
3 

fe = e,/i,T 



pAL=2 0fc) pAL=0(ifc) 
^ Yukawa ^ ' 



a=l 



Yukawa 



piy (ak) pW (ak) 



A^a \ Yukawa 



pAL=2(fcj) _ pAL=0(A;i) 
^ Yukawa ^ ^ 



Yukawa 



pD (aO 



n. 



21 77AL, 



(3.16) 

(3.17) 
(3.18) 
(3.19) 



VZ 21 riALi 



V 63 r^ALj 9 TfZ VAL, J 



(3.20) 



'm. ,4 ?7AL \ N^L , 1^1 , 1 ^AL 5 77Ar„ 77AL A N^uO 

.vZ ^ 21 valJ^^^' + I' + 9 ,A., + 63 ,A.J ^^'^^ 



, A , 5 77AL 1 VN^ VAL \ N^Q 

+ 1 + 7^ ^ n ^Lu 

V 63 r]ALi 9 Tyjv}'^ TyAL,/ 



(3.21) 
(3.22) 



1 



vZ 21 r^AL, 



1 + 777 7l$ 

21 t^alJ 



+ 1,' 21^;^;^-^'^ 



(3.23) 



4 77Al\ /v.L 



4 77 AL \ Af^^Vu 



21 ?7AL 

' 4 riAL \ N^^t 

.^21 ^^„^alJ^^^'^ 



21 rjALi 



7l$ 



(3.24) 
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63 r)ALi 



(3.25) 



(al) 
Gauge 



1 



2 + 



21 r/AL,, 



+ 1 + :^ 



Val 



N L 



+ 1 + 



21 ^AL, 



7lI 



21 r;AL< 



7lk 



-, AL=2 (jfe) 
Yukawa 



AL=0 (jfc) 
Yukawa 



1 

1 



1 + 



1 + 



21 r^AL,- 



^?AL 



21 r/AL, 



/Lj<S> Lj<S>t LjL^ 

7Lfc$ + Ti.^t + 7$$t 



(3.26) 



(3.27) 



(3.28) 



Notice that the would-be singularities in the hmit of a vanishing t^al; in (3.18)-(3.21) and 
(3.23)-(3.28) are exactly cancelled by corresponding factors r]ALi that multiply the collision 
terms in the BEs (3.15) and (3.16). We should also note that the flavour-diagonal AL = 
processes, with k — j, do not contribute to the BEs, as it can be explicitly checked in (3.16). 
Finally, it is worth commenting on the earlier form of the BEs, obtained in [17]. This can 
be recovered from (3.15)-(3.28), after summing over the three lepton-doublet BEs, with 
the assumption that riALj = I^t-al, and after neglecting SM chemical potential corrections. 



4 Out of Equilibrium Sphaleron Effects 

In the SM, the combination of the baryon and lepton numbers, B + L, is anomalous [49]. 
Although at low energies this B + L violation is unobservably small, at temperatures close 
to and above the electroweak phase transition, e.g. for T > 150 GeV, thermal fluctuations 
more efficiently overcome the so-called sphaleron barrier allowing rapid B + L violation in 
the SM [5]. 

The temperature dependence of the rate of S -|- L violation is of particular interest 
in models of low-scale leptogenesis. Any lepton asymmetry produced after the {B -\- L)- 
violating interactions drop out of thermal equilibrium will not be converted into a baryon 
asymmetry. Therefore, in electroweak-scale leptogenesis scenarios, it is important to con- 
sider the rate ol B -\- L violation in the BEs, in order to offer a more reliable estimate of 
the final baryon asymmetry. 

The rate of {B -\- L)-violating transitions has been studied in detail in [33, 50] for 
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temperatures satisfying the double inequality 

Mw{T) < T < ^^^wiJl 

where = 5'^/47r is the 811(2)^ fine structure constant, Mw{T) — gv{T)/2 is the T- 
dependent VF-boson mass and 

v{T) = ^(0)(^l-^y (4.2) 

is the T-dependent VEV of the Higgs field. The critical temperature of the electroweak 
phase transition, T^, is given by [51] 



1 , ag' , , h 



where A is the quartic Higgs self-coupling, g' is the U(l)y gauge coupling and ht is the 
top-quark Yukawa coupling. 

The rate of S -|- L violation per unit volume is [50] 



r 00. 



IHB^L) ^ ^ = ^N^ANVUi^] afe-^^^'^n. (4.4) 



V 2 71 4 



According to (4.1), this expression is valid for temperatures T < Tc. The various quantities 
in (4.4) are related to the sphalcron dynamics and are discussed in [33,50]. Following the 
notation of [50], the parameters a;_, A/tr and AArot are functions of X/g^, Kot = Stt^ and 
0:3 = g'i/4:'K, where 

Egp is the energy of the sphaleron and is given by 

E^^a'-^, (4.6) 

aw 

where A is a function of X/ g'^ and is of order 1 for all phenomenologically relevant values 
oi X/g"^. The dependence of the parameter k on X/ g"^ has been calculated in [33,50], using 
various techniques. The results of those studies are summarized in Table 1, where the 
values of k and the other sphaleron-related parameters are exhibited for X/g"^ — 0.278, 
which corresponds to a SM Higgs-boson mass Mh of 120 GeV. 

Given the present experimental hmits on the SM Higgs-boson mass, Mh ^115 GeV, 
it can be shown that the electroweak phase transition in the SM will either be a weakly first 



19 









Mr 




A 


0.278 


0.806 (gv) 


11.2 


7.6 


0.135 - 1.65 


1 



Table 1: Values of the various parameters in (4.4) for X/g"^ — 0.278, corresponding to a 
SM Higgs-boson mass of 120 GeV. 

order one, or even a second or higher order one, without bubble nucleation and the forma- 
tion of large spatial inhomogeneities in particle densities. Therefore, the use of a formalism 
describing the (B + L) -violating sphaleron dynamics in terms of spatially independent B- 
and L- number densities ub and riLj niay be justified. Further refinements to this approach 
will be presented elsewhere. 

We should bear in mind that heavy Majorana neutrino decays, sphaleron effects and 
other processes considered in the BEs (3.15) and (3.16) act directly on the number densities 
of SU(2)i lepton doublets, nAL^^ .^. However, the quantity usually referred to as lepton 
number, L, has a contribution from the right-handed charged leptons le,n,TR as well. In 
thermal equilibrium, one may relate the asymmetry in right-handed charged leptons to the 
asymmetry in SU(2)i lepton doublets by virtue of (3.9), leading to the result 

?7A/« = 2 ^^^i ~ 21 ^^-^ ■ ^ 
The change in lepton flavour can be thought of as having two components, one component 
termed leptogenesis due to lepton-number- violating processes considered in Section 3, and 
another due to the [B + L)- violating sphalerons: 



dz dz 
where, up to 0{v/T) corrections. 



+ 

leptogenesis 



dr]i 



dz 



(4.8) 

sphaleron 



dz 



3 dr]^Li 2 d77AL 



, , . 2 dz 21 dz 

leptogenesis 



The BEs determining the leptogenesis component of (4.8) have been discussed in Section 3. 
We shall now discuss the BEs determining the sphaleron component of (4.8), and the 
generation of a net S-number asymmetry. 

Within the context of the above formalism, the sphaleron components of the BEs for 
Ub and riL^ are given by 

dt ^ 
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^^^3 _^ 2>HnL = ( el3{i^BQB+^LQL} _ ^I3{ij.bQ'b+i^lQ'l} 
dt ^ ^ 

_ ^-P{nBQB+^^LQL) + ^-PMs+i^LQ'-^) ^ 7^^^^^) , (4.IO) 

where uq is the number of generations and (5 — 1/T. Furthermore, Qb{l) is the baryonic (or 
leptonic) charge of the system before the (S+L)-violating sphaleron transition and Q'b{l) 
the charge after the transition. Khnkhamer and Manton showed [52] that a sphaleron carries 
a baryon (and lepton) number of 71(3/2, therefore Qb — Q'b — t^g/'^ and — Q'hi = 1/2. 
Finally, assuming that the baryon and lepton chemical potentials are small with respect to 
the temperature, the BEs (4.10) may be approximated by 



dns 
"df 

dni^ 
dt 



i 



Notice that the BEs (4.11) are linear in the chemical potentials, which is a very useful 
approximation for our numerical estimates. 

We now need to determine the relation between the baryon and lepton chemical po- 
tentials and their respective number densities. These relations can be found by considering 
the effective potential, V , of the Higgs and the 811(2)^ and U(l)y gauge fields. They have 
been computed in [30] at finite temperatures, for small chemical potentials, //s, jiL <^ T 
and when v{T) < (a few) x T. In this framework, the neutrality of the system with re- 
spect to gauge charges can be accounted for by minimizing the potential with respect to 
the temporal components of the SU(2)i and U(l)y gauge fields, Wq (a = 1,2,3) and Bq, 
respectively. The baryon and lepton number densities are then given by 

dV dV ,^ , 

= , = ■ 4.12 

For the SM with 3 generations and 1 Higgs doublet, we obtain 

77T^ + 27v\T) ^ 22T^ + 3v^{T) ^ 



51^2 



fsiriL, - 3nB + 4: ^ ulA 

j =e,^J■,T 

(^SuB - ^Y^^Lj) ■ (4-13) 



J =e,^J■,T 

484 



153(4472 + 17 v2(r 



J =e,fi,T 
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Employing the relations (4.13), we may now extend the system of BEs (3.15) and 
(3.16), by explicitly taking the {B + L)-violating sphaleron transitions into account. 



dr]B 
dz 



H(z = 1) 



(4.14) 



X 



28 ^ 225 v^m 

j =e,iJ,,T 



561 T2 



108 ^ 

3 =e,IJ.,T 



A(B+L) 



drji 



drji 



dz 



dz 



leptogenesis 



+ 1 
3 dz 



with 



(3366/7r2) T2 -/^(^b+l) 



(4.15) 



(4.16) 



A{B+L) T2 + 51^;2(r) 

The leptogenesis component of (4.15) may be determined using relation (4.9), along with 
the BEs (3.15) and (3.16). 

Observe that in the limit of infinite sphaleron transition rate, V /^(b+l)/ H{z = 1) ^ 
cxo, and at high temperatures T ^ v{T), the conversion of lepton-to-baryon number den- 
sities is given by the known relation: 



28 ^ 

3 =e,H,r 



(4.17) 



To account for the T-dependent (5+L)- violating sphaleron effects, our numerical estimates 
given in the next section will be based on the BEs (3.15), (3.16), (4.9), (4.14) and (4.15). 



5 Numerical Examples 

We shall now analyze RrL models that comply with the constraints obtained from the 
existing low-energy neutrino data [11, 53] and provide successful baryogenesis. As was 
discussed in Section 2, our specific choice of model parameters will be motivated by the 
naturalness of the light and heavy neutrino sectors. 

Phenomenologically relevant RrL models can be constructed for an SO (3) invariant 
heavy neutrino mass of the size of the electroweak scale, e.g. niN — 250 GeV [cf. (2.3)], if 
\a\ ~ |6| |c| and |a|, \b\ ~ 10~^. To protect the r-lepton number from wash-out effects, 
we also require that the small U(l);-breaking parameters |£^e,/i,T| be no larger than about 
10~^ and |c| < 10~^. For definiteness, the model parameters determining the hght neutrino 
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a/(mjv)5 (GeV)-^ 


(AM5)22/mjv 


(AM5)33/miv 


(AM5)23/miv 


6.0 X 10-^ 


4.0 X lO-'' 


5.2 X lO-'^ 


(6.8- 0.6 i) X 10-9 



Table 2: Choices of the parameters a, {AMs)22, (AM5)33 and {AMs)23, consistent with 
hght neutrino data. 

mass spectrum are chosen to be (in arbitrary complex units) 

AniN 2 . n 21 13 49 ,^ 

= 4, se ^ 2 + = 77^ > = ' 

rnisf 250 50 128 

where the ratio b/a is kept fixed: 



(5.2) 



6 _ 19 
a ~ 50 ■ 

The actual values selected for the relevant parameters a, (AM5)22, (AM5)33 and (AM5)23 
vary with the S0(3) invariant mass mjv. As we will see in more detail below, Table 2 
illustrates choices of these parameters consistent with the light neutrino data. For ttin 
in the range 100-1000 GeV, the chosen parameters are consistent with the naturalness 
condition mentioned in Section 2, whilst giving rise to phenomenologically rich models. 

In our numerical analysis, we will focus on 4 examples, with m^r = 100, 250, 500, 
and 1000 GeV. Clearly, the model parameters selected in (5.1), (5.2) and Table 2 imply 
that all the scenarios have the same tree-level hght neutrino mass matrix: 

/-1.27 3.63 2.96 \ 

3.63 1.89 0.370 
y 2.96 0.370 -0.544 y 

This leads to an inverted hierarchy of light neutrino masses, with mass differences and 
mixings compatible with the current 3a bounds [53]. Adopting the convention rriiy^ < 
rUy^ < m^2, we find the mass squared differences and mixing angles 

ml^ - m\ = 7.54 x 10"^ eV^ , m\ - ml^ = 2.45 x 10"^ eV^ , 

sin2^i2 = 0.362, sin^ ^23 = 0.341, sin^ ^13 = 0.047. (5.4) 

Since the mass matrix (2.8) is rank 2, one light neutrino will be massless at the tree level 
(m,^3 = 0) , thus fixing the absolute scale of the light neutrino hierarchy. 

The remaining soft S0(3)-breaking parameters, {AMs)n, {AMs)i2, (AM5)i3, do not 
affect the light neutrino mass spectrum. These together with the parameter c play a key 



X 10-^ eV. (5.3) 
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role in obtaining the correct BAU and are exhibited in Table 3. We choose {AMs)ii to be 
relatively large, {AMs)ii ~ lO^^mjv, providing large mass differences ImN^—mN-^l/mN and 



|mjV3 — mTvJ/mTv ~ fO^"''. Such a choice is consistent with thermal and RG effects running 
from the GUT scale ~ fO"*^^ GeV to the electroweak scale ~ mjv (see also the discussion 
in Section 2). The other two soft S0(3)-breaking parameters, {AMs)i2 and {AMs)i3, are 
selected so as to give the observed BAU. 

To assess the degree of cancellation between tree-level and RG contributions to AMs, 
we introduce the parameter rc defined as 

In (5.5), the product is taken over contributions where |(AM^^)jj| > \{AMs)ij\. 

The parameter tq is always greater than 1 and represents that the degree of cancellation is 
1 part in rc- From the values of rc displayed in Table 3, we observe that electroweak-scale 
heavy Majorana neutrinos are favoured by naturalness. 

The baryon asymmetry predicted for each model can be determined by solving the 
EEs (3.15), (3.16), (4.9), (4.14) and (4.15), and using the collision terms derived in Ap- 
pendix A and [17]. These coUision terms are calculated in the basis where the charged- 
lepton and heavy-Majorana mass matrices are positive and diagonal. They have been 
appropriately expressed in terms of the one-loop resummed effective Yukawa couplings 
derived in [17]. It is worth noting that all SM reactions, including those involving the 
e- Yukawa coupling, are in full thermal equilibrium for the temperatures relevant to our 
scenarios, T < 10 TeV [51,54]. Moreover, since heavy Majorana neutrino decays are ther- 
mally blocked at temperatures T > 3mjv„ [55], we will only display numerical estimates of 
the evolution of lepton and baryon asymmetries, for z — rriN^/T > 0.1. Nevertheless, as 
we will see below, the predictions for the final BAU are relatively robust in RL models, 
because of the near or complete independence on the primordial baryon and lepton number 
abundances. 

Some of the Yukawa and gauge-mediated collision terms contain IR divergences, which 
are usually regulated in thermal field theory by considering the thermal masses of the 
exchanged particles [43]. To assess the theoretical errors introduced by the choice of a 
universal thermal mass regulator (see the discussion in Appendix A), we have estimated 
the response of the final baryon asymmetry under variations of the IR mass regulator miR. 
We find that the predicted BAU only varies by ±7 %, for a variation of miR by ±25 %. 
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Figure 2: The predicted evolution of rj^^ and tjb, for models with m^v = 100, 250, 500 
and 1000 GeV, and rj'^^ = 1. The model parameters are given in (5.1), (5.2), and Tables 2 
and 3. The horizontal grey dashed line shows the baryon asymmetry needed to agree with 
observational data. 
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mjv (GeV) 


(AM5)n/mjv 


(AM5)i2/miv 


{^Ms)n/mN 


c 




100 


1.0 X 10-^ 


-1.00 X 10-*^ 


-5.5 X 10-^° 


1.0 X 10- 


-7 


5 


250 


1.0 X 10-^ 


-1.36 X 10-*^ 


-8.0 X 10-^° 


1.5 X 10- 


-7 


39 


500 


1.0 X 10-^ 


-1.36 X lO-'^ 


-8.8 X 10-^° 


2.0 X 10- 


-7 


264 


1000 


1.0 X 10-^ 


-1.9 X lO-'' 


-1.0 X lO-'^ 


2.5 X 10- 


-7 


1240 



Table 3: Choices of {AMs)u, {AMs)i2, {AMs)i3 and c, which, in conjunction with those 
in Table 2, lead to successful baryogenesis. rc parameterizes the degree of cancellation 
between radiatively induced and tree-level contributions to AM5. 

The BEs are solved numerically, using the Fortran code LeptoGen^ Fig. 2 shows 
the predicted evolution of the baryon and individual lepton asymmetries, tjb and rj^, 
as functions of the T- related parameter z — uini/T, for each of the 4 examples, with 
tun — 100, 250, 500 and 1000 GeV. The specific model parameters are given in (5.1), 
(5.2), and Tables 2 and 3. Each scenario had an initially thermal heavy Major ana neutrino 
abundance and zero initial baryon and lepton asymmetries, i.e. 77™^ = 1 and 77^ = rj^ — 0. 
The 4 panels show that the large asymmetry is slightly reduced by less significant, but 
opposite sign Lg and Lf^ asymmetries. Clearly visible in each scenario is the effect of the 
rapidly decreasing rate oi B + L violation; the lepton and baryon asymmetries quickly 
decouple at T ~ T^. This decoupling is particularly pronounced in the uin — 100 GeV 
scenario where the baryon asymmetry freezes out exactly when the lepton asymmetry is 
maximal. In particular, the rapid decoupling of tjb from r)Li at temperatures T close to 
Tc has the virtue that, unlike rjii, rjB remains almost unaffected from ordinary SM mass 
effects due to a non-zero VEV v{T) [cf. (4.2)], since it is v{T ~ < v{T = 0). 

Fig. 3 shows the evolution of the baryon asymmetry for varying initial lepton, baryon 
and heavy neutrino abundances. For the 250 GeV scenario. Fig. 3(a) illustrates the near 
independence of the resultant baryon asymmetry on the initial conditions. Even for the 
most extreme initial conditions r]^ = =f0.1 and 77^ = ±0.1, the variation in the final 
baryon asymmetry is only ± 38%. 

For heavy neutrino masses m^y < 250 GeV, the dependence on initial conditions be- 
comes stronger. In the m^r = 100 GeV scenario, Fig. 3(b) shows the dependence of the final 
BAU on the initial lepton and baryon asymmetries in a RrL scenario with rriN = 100 GeV. 
It is interesting to observe that the final B asymmetry will remain almost unaffected, even 
tLeptoGen may be obtained from http://hep.man.ac.Uk/u/thomasu/leptogen 
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Figure 3: The (in) dependence of the final baryon asymmetry on the initial lepton, baryon 
and heavy Majorana neutrino abundances for (a) = 250 GeV and (b) = 100 GeV. 
The model parameters and the meaning of the horizontal grey line are the same as in Fig. 2. 
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T 

Figure 4: The complete independence of the final BAU on the initial lepton and baryon 
abundances for the rriN = 500 GeV scenario, with model parameters the same as in Fig. 2. 

if tlie primordial baryon asymmetry 77^ at T > 10 m^v is as large as 10"^, namely two orders 
of magnitude larger than the one required to agree with observational data. 

In RrL scenarios with ttin > 250 GeV, the final baryon asymmetry is completely 
independent of the initial conditions. This is illustrated in Fig. 4 for the RrL scenario with 
m^v = 500 GeV. In this numerical example, it is most striking to notice that the prediction 
for the final BAU remains unchanged, even if the initial conditions are set at temperatures 
helow the heavy neutrino mass scale m^v, e.g. at T ~ 0.5 mAr. 

Some insight into the independence on initial conditions is provided by Fig. 5. The 
ratios of various collision terms to the Hubble parameter are plotted for the rriN = 250 GeV 
scenario. These ratios show that RL can take place almost completely in thermal equilib- 
rium; in certain cases, the reaction rates are many orders of magnitude above the Hubble 
parameter H{z =1). In spite of this fact, RL (RrL) can successfully generate the required 
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Figure 5: The ratio of various collision terms to the Hubble parameter plotted for the 
TTiAT = 250 GeV scenario presented in Fig. 2. 

excess in L (Lt), because of tlie resonantly enlianced CP asymmetry. 

To allow for a simple quantitative understanding of the baryon asymmetry in RrL 
(and similar) scenarios, we need to introduce the individual lepton flavour i^-factors 

^iv„ - Hiz = l) ■ ^^-^^ 

Note that the decay widths are calculated in terms of the one-loop resummed effective 
Yukawa couplings [17]. 

Table 4 shows the various components of K]^^ for the = 250 GeV scenario. 
This explicitly demonstrates how the texture provided by (2.5) and (2.7) allows for a 
heavy Majorana neutrino to decay relatively out of equilibrium, whilst simultaneously 
protecting the r- lepton number from being washed-out, even though large e- and //-Yukawa 
couplings to A''i_2 exist. Bear in mind that we use the convention ttini < < upon 
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Table 4: Individual lepton flavour K -factors for the uin — 250 GeV scenario. 



diagonalization of the heavy Majorana neutrino mass matrix Ms- As can be seen from 
r^^'^ ~ 10-100 and a CP-asymmetry SJ^^ 



Table 4, i^-factors K'^'^''^ ~ 10-100 and a CP-asymmetry ~ ~10 ^ are sufficient to 



generate a large r-lepton asymmetry. Although the fT-factors K'f^'^,^ associated with 
and the e and // leptons are enormous of order lO'^-lO^'^, these turn out to be harmless to the 
T-lepton asymmetry, as the latter is protected by the low r-lepton i^T-factors KJ^^^^a ^ 

An order of magnitude estimate of the final baryon asymmetry, including single lepton 
fiavour effects, may be obtained using 

3 Ty-l 

rjB 10-2 X E E e-("*^«-"^~i)/'"^i 4 iT^r- ■ (5.7) 

1=1 " 

The above estimate for rjB is also consistent with the one stated earlier in [26]. In (5.7), 
the i^-factors are summed in the following way: 

Kn^ - E K^n^, Ki ^Y. e-(™-"--i)/'"-i K'j,^ . (5.8) 

1=1 No, 

Notice that all i^'-factors are evaluated at T = m^r-^ (i.e. z = tiin^/T = 1), where mAr-^ is 
the lightest of the heavy Majorana neutrinos. The intuitive estimate (5.7) is applicable for 
all leptogenesis scenarios satisfying the approximate inequality 

KiNo^ ^ 1, (5.9) 

for each of the lepton flavours / and the heavy Majorana neutrinos A^i,2,3- The inequal- 
ity (5.9) ensures that the energy scale m^Vi can be identified as the true scale of leptogenesis. 

In RL scenarios, such as RrL, the importance of taking individual lepton fiavour 
effects into account can be demonstrated by comparing (5.7) with the naive estimate, in 
which lepton fiavour effects are treated indiscriminately in a universal manner, 

rjl"""- - -10-2 X e-^"'^"-"'^!)/"*^! — , (5.10) 

N ^ 
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Figure 6: The predicted evolution of ?7Ltotai Vn^ ^or a model with a heavy neutrino 
spectrum: m^N^ = 3mArj, m^^ = 2m7Vi and m^^ = lO"*^" GeV. 

where K = J2i=e.fiT ^i- the RrL scenario with m^r = 250 GeV, the dominant contribu- 
tion to this estimate will come from N^, with a total CP asymmetry 6^.^ ~ 10~^. Taking 
the ratio of the two estimates yields 



„ univ. 

Vb 


6^3 






|c|2 


Vb 




KkK 




|aP + |6|2 



Thus, without considering single lepton flavour effects in this particular RrL model, one 
obtains an erroneous prediction for the BAU, which is suppressed by 6 orders of magnitude 
and has the wrong sign. These estimates are confirmed by solving the total lepton number 
BEs presented in [17]. 

In a hierarchical scenario, the number densities of the heavier neutrinos A^2,3 at T = 
mTVi will be Boltzmann suppressed. To account for this phenomenon, we have included the 
Boltzmann factors e"^™^""™'^!^/'"^! in the estimates (5.7), (5.10) and in the definition of 
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Ki. Clearly, in RL models with each heavy neutrino nearly degenerate in mass, this last 
factor can be set to 1. 

Flavour effects can also play a significant role in mildly hierarchical scenarios. Figure 6 
shows the predicted evolution of the lepton asymmetry in a scenario where mjvg = Sm^r^, 
171^2 = 2m7Vi "^JVi = 10"^° GeV. The Yukawa texture was chosen to be consistent with 
light neutrino data and a normal hierarchical light neutrino spectrum. In this example, 
neglecting individual lepton flavour effects introduces an 0(10) suppression of the flnal 
lepton and baryon asymmetry. 

In fully hierarchical scenarios satisfying (5.9), it can be seen that the estimates (5.7) 
and (5.10) are completely equivalent. A large hierarchy in the heavy neutrino spectrum, 
combined with the condition (5.9), implies that the final lepton asymmetry is determined 
entirely by the decay of the lightest heavy Majorana neutrino A^i. This fact makes it 
impossible for a single lepton fiavour to be protected from wash-out, whilst the neutrino 
decays out of equilibrium. 

Likewise, in fiavour universal scenarios, where i]Le,^,T = \vl, the estimates (5.7) and 
(5.10) are completely equivalent for both nearly degenerate and hierarchical leptogenesis 
scenarios. 

Our numerical analysis presented in this section has explicitly demonstrated that 
models of RrL can provide a viable explanation for the observed BAU, in accordance with 
the current light neutrino data. In the next section, we will see how the scenarios considered 
here have far reaching phenomenological implications for low-energy observables of lepton 
fiavour/number violation and collider experiments. 

6 Phenomenological Implications 

RL models, and especially RrL models, can give rise to a number of phenomenologically 
testable signatures. In particular, we will analyze the generic predictions of RrL models 
for the OuPP decay, and for the LFV processes: fi —>■ ej, fi —>■ eee and fi —>■ e conversion 
in nuclei. Finally, we will present simple and realistic numerical estimates of production 
cross sections of heavy Majorana neutrinos at future e^e~ and /i^/i^ colliders, and apply 
these results to the RrL models. 
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6.1 Ou(3(3 Decay 



Neutrinoless double beta decay {Qv(5(5) corresponds to a process in which two single (5 
decays [56-58] occur simultaneously in one nucleus. As a consequence of this, a nucleus 
{Z, A) gets converted into a nucleus {Z + 2, ^4), i.e. 

zX ^ Z+2X + 2e . 

Evidently, this process violates L-numbcr by two units and can naturally take place in 
minimal RL models, in which the observed light neutrinos arc Majorana particles. The 
observation of such a process would provide further information on the structure of the 
light neutrino mass matrix m". 

To a very good approximation, the half life for a Oi/(3(3 decay mediated by light 
Majorana neutrinos is given by 

[^1/?']"' - ^|A^o.,,rGoi, (6.1) 

where (m) denotes the effective Majorana neutrino mass, mg is the electron mass and 
M Qyjip and Goi denote the appropriate nuclear matrix element and the phase space factor, 
respectively. More details regarding the calculation of Ty2^ may be found in [56-58] . 

In models of interest to us, the effective neutrino mass can be related to the entry 
{11} (= {ee}) of the light neutrino mass matrix m'' in (2.8), i.e. 



|(m)| = l(m^), 



,2 



Am AT o n 

a - e t. 

rriN 



(6.2) 



2mjv 

As has been explicitly demonstrated in the previous section, RrL models realize a light 
neutrino mass spectrum with an inverted hierarchy [59], thus giving rise to a sizeable 
effective neutrino mass. The prediction for | (m) | in these models is 

|(m)| = |(m'^)ee| ~ 0.013 eV. (6.3) 

Such a prediction lies at the very low end of the value of the effective Majorana neutrino 
mass, reported recently by the Heidelberg-Moscow collaboration [60]. There are proposals 
for future 0z//9/9-decay experiments that will be sensitive to values of | (m) | of order 10~^ [36], 
significantly increasing the constraints on this parameter. 
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Figure 7: Feynman graphs pertaining to the effective 'ye/j, and Ze/j, couphngs. 
6.2 ij, ^ ej 

As shown in Fig. 7, heavy Majorana neutrinos may induce LFV couphngs to the photon 
(7) and the Z boson via loop effects. These couplings give rise to LFV decays, such as 
II ^ ej [61] and // — > eee [62]. Our discussion and notation closely follows the extensive 
studies [62,63]. Related phenomenological analyses of LFV effects in the SM with singlet 
neutrinos may be found in [64-67]. 

To properly describe LFV in low-energy observables, we first introduce the so-called 
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Langacker-London (LL) parameters [68]: 

= 1- E ^ {rn*^M*s-'Ms'ml)^^, (6.4) 

i'=e,//,T 

where ttid — h"'^ v and Bi^i^ are mixing-matrix factors close to 1 that multiply the SM 
tree-level Wlvi' vertices [39]. The LL parameters (s^')^ quantify the deviation of the actual 
squared Wlui' couplings (summed over all light neutrinos) from the corresponding sum 
of squared couphngs in the SM. The parameters (s^'''''^)^ are constrained by LEP and 
low-energy electroweak observables [68,69]. Independent constraints on these parameters 
typically give: (s^'''^'^)^ ^ 10~^. As we will see in a moment, however, LFV observables 
impose much more severe constraints on products of the LL parameters, and especially 
on sl^sl^. 

As a first LFV observable, we consider the decay — > 67. The branching fraction for 
the decay — > 67 is given by 

where = 2.997 x 10~^^ GeV [11] is the experimentally measured muon decay width, and 
is a composite form-factor defined in [62]. In arriving at the last equality in (6.5), we 
have assumed that 171% ^ M^r, for a model with two nearly degenerate heavy Majorana 
neutrinos. In this case, one finds that GUj^ ^ "SiTsl , where 4> is an unobservable model- 
dependent phase. Confronting the theoretical prediction (6.5) with the experimental upper 
hmit [11] 

5exp(/i^e7) < 1.2 X 10-'\ (6.6) 
we obtain the following constraint: 

sl^sl^ < 1.2 X 10-^ (6.7) 

This last constraint is stronger by one to two orders of magnitude with respect to those 
derived on {s^j^y and (s^'')^ individually. 

In RtL models, only two of the right-handed neutrinos, i'2r and u^r, which have 
appreciable e- and /x- Yukawa couplings, will be relevant to LFV effects. In this case, the 
LL parameters {s'^)^ and {s'l^Y are, to a very good approximation, given by 

K-r - ^ , isin' - ^ ■ (6.8) 
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Then, the following theoretical prediction is obtained: 

Bit, ^ 67) = 9 • 10-^ X . (6.9) 

For the particular scenarios considered in Section 5, we find — > 67) ~ 10~^^. These 

values are well within reach of the MEG collaboration, which will be sensitive to — > 
67) ~ 10-1^ [34]. 



6.3 /J, eee 

As illustrated in Fig. 8, quantum effects mediated by heavy Majorana neutrinos may also 
give rise to the 3-body LFV decay mode — > e~e'^e~. The branching ratio for this LFV 
decay may conveniently be expressed as 



B{fj, ^ eee) — 



OL. 



mj^ m 



24 5 767r3 



+ 44|Fr - |2 + 164 Re[(Fr + XlGit' 



■484 Re (Ff - F^^^)Gf * + 324|G^;^f ( In-f - — j ^ . (6.10) 



The expressions F!^, F^^ and F^^^ are composite form-factors, defined and computed 
in [62]. In the limit » and up to an overall physically irrelevant phase factor e**^, 
these composite form-factors simphfy to [62] 
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(6.11) 

(6.12) 
(6.13) 



Correspondingly, the analytic result (6.10) in the same limit may be cast into the form: 



S(/i — >• eee) ~ 
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(6.14) 



It can be seen from (6.14) that the so-called non-decoupling terms proportional to m%/M^ 
are always multiplied with higher powers of the LL parameters. In general, these terms 
do not decouple and become very significant [62], for large heavy neutrino masses niN and 
fixed values of s^' , which amounts to scenarios with large neutrino Yukawa couplings | h^f \ ^ 
0.5 [39]. However, these non-decoupling terms are negligible, as long as m^/Mw <C 1. 
This is actually the case for the RrL models discussed in Section 5. Neglecting terms 
proportional to m%/My^ and m^/M^, we may relate B{ii ^ eee) to B{fi ^ ej) through: 



B{ii eee) ~ 8.2 • 10"^ x 



0.8 In 



m 



N 



+ 0.5 In' 



m 



N 



B{fi ^ 67) . (6.15) 



For example, for an RrL model with ttin — 250 GeV, (6.15) imphes 



Bit, 



eee] 



1.4 • 10-2 X B{ti 67) ~ 1.4 • 10 



-14 



(6.16) 



This value is a factor ~ 70 below the present experimental bound [11]: -Bexpl/^ — ^ eee) < 
1.0 X 10~^2. In this respect, it would be very encouraging, if higher sensitivity experiments 
could be designed to probe this observable. 



6.4 Coherent ji — > e Conversion in Nuclei 

One of the most sensitive experiments to LFV is the coherent conversion of /i ^ e in nuclei, 
e.g. ji~ '^Ti — > e~ ^Ti [70,71]. The Feynman graphs responsible for such a process are 
displayed in Fig. 8. 

Our calculation of /x — e conversion in nuclei closely follows [63,70,71]. We consider 
the kinematic approximations: ~ ~f^^\ and Pg ~ |pe| ~ T^in which are valid for /i — > e 
conversion. Given the above approximation, the // — > e conversion rate in a nucleus with 
nucleon numbers ( A^, Z) , is given by 

o rAAZl - r[/x(iV,Z)^e(iV,Z)] ^ alr^ajml Z^^ 

^-^"^'^^ = r[MA^,^)- capture] ~ 32.^M^r..p, ^ "^^-"^'^^l ' ^'-^'^ 
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Figure 8: Feynman graphs responsible for /i — > e e"^e (// — > e conversion in nuclei). 

where ctem = 1/137 is the electromagnetic fine structure constant, Z^s is the effective 
atomic number of coherence and Fcapt. is the muon nuclear capture rate. For 22Ti, ex- 
perimental measurements give Zes ~ 17.6 for ||Ti [72] and r[// I^Ti — > capture] 
1.705 X 10-^^ GeV [73]. Moreover, \F{-ml) \ ^ 0.54 is the nuclear form factor [74]. Finally, 
Qw — Vu{2Z + N) + Vd{Z + 2N) is the coherent charge of the nucleus, which is associated 
with the vector current. Its explicit form is given by 

K = ^4 - G^,^ - Ff) + i (Ff - FC-) , (6.18) 

= - ^ 4 (^r - - ) - J (i^r + ^^b^:'') ■ (e.w) 

The composite form-factors F^^^^ and F^^^'^ are defined in [63]. In the SM with two nearly 
degenerate heavy Major ana neutrinos and in the limit m%/M^r 3> 1, these form- factors 
can be written down in the simplified forms: 

In the same limit m%/M^^ S> 1, B^e{N, Z) is given by 
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rriN (GeV) 


B{fi 67) 


B{fi eee) 


5^e(26,22) 


100 


6.2 X 10-12 


3.8 X lO-i'i 


9.2 X 10-13 


250 


9.9 X 10^13 


1.4 X 10-1^ 


4.5 X 10-13 


500 


2.5 X 10-13 


9.7 X 10-15 


2.0 X 10-13 


1000 


6.2 X 10-1^ 


4.9 X 10-15 


7.7 X 10-1^ 



Table 5 : Branching fractions for tlie 3 LFV processes ji — > 67, ji — > eee and coherent 
II ^ e conversion in g^Ti nuclei. 
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For the 22Ti case, S^e(26, 22) is related to S(/x — > 67) through 



S^e(26,22) ~ 0.1 X 



1 + 0.5 In 



ra 



N 



B{ii^ 67) 



(6.22) 



On the experimental side, the strongest upper bound on S^e(A^, Z) is obtained from 
experimental data on // — > e conversion in f|Ti [75] : 



(26,22) < 4.3 X 10-1% 



(6.23) 



at the 90% CL. However, the proposed experiment by the MECO collaboration [35] will be 
sensitive to conversion rates of order 10- 1^. 

In the RtL model with ttin — 250 GeV, one obtains, on the basis of (6.22), the 
prediction for // — > e conversion in ^2^1: 



5^e(26,22) ~ 0.46 X B(/x ^ 67) ~ 4.5x10" 



-13 



(6.24) 



The above prediction falls well within reach of the sensitivity proposed by the MECO 
collaboration. 

In Table 5, we summarize our results for the branching ratios of the 3 LFV processes: 
II — > 67, n — > 666 and coherent 11 ^ e conversion in 22Ti nuclei, for each RrL model 
considered in Section 5. 
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Figure 9: Diagrams related to the production of the heavy Majorana neutrinos A^2,3 at 
future e'^ e~ {ii^ high-energy colliders. 

As a final general remark, we should mention that RrL models, and leptogenesis 
models in general, do not suffer from too large contributions to the electron electric dipole 
moment (EDM) [12,76], which first arises at two loops. The reason is that EDM effects 
are suppressed either by higher powers of small Yukawa couplings of order 10^"^ and less, 
or by very small factors, such as {itlm^ — mN2z)/f^N ~ 10^'^. The latter is the case in RrL 
models, which leads to unobservably small EDM effects of order 10"^'^ e cm, namely 10 
orders of magnitude smaller than the present experimental limits [1 1] . 



6.5 Collider Heavy Majorana Neutrino Production 

If heavy Majorana neutrinos have electroweak-scale masses and appreciable couplings to 
electrons and muons they can be copiously produced at future e+e~ [77, 78] and 
colliders. As shown in Fig. 9, this is exactly the kinematic situation for the heavy Majorana 
neutrinos iV2,3 described by the RrL models. The heavy Majorana neutrino A^i has a very 
small coupling to leptons and it would be very difficult to produce this state directly. 

For collider c.m.s. energies ^/s ^ ttt-at, the t-channel VF~-boson exchange graphs 
will dominate over the Z-boson exchange graph, which is s-channel propagator suppressed 
(see Fig. 9). In this high-energy limit, the production cross section for heavy Majorana 
neutrinos approaches a constant [77], i.e. 

a[e+e-(/.V)-A^2,3^] ^ ^('L'^'f ^ lOfhx (^^j . (6.25) 

Since s'^ in RrL models, the produced heavy Majorana neutrinos A^2,3 will have the 
characteristic signature that they will predominantly decay into electrons and muons, but 
not into r leptons. Assuming that uin > Mh, the branching fraction of A^2,3 decays into 
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charged leptons and into bosons decaying hadronically is 

s(7V2,3^e^,//^W^±Hjets)) 1 x ^ = ^ . (6.26) 

Given (6.25), (6.26) and an integrated luminosity of 100 fb~^, we expect to be able to 
analyze about 100 signal events for {s^^''^)'^ — 10~^ and ijin ~ 300 GeV, at future e+e~ 
and /jl'^/jl' colliders with c.m.s. energy ^/s — 0.5-1 TeV. 

These simple estimates are supported by a recent analysis, where competitive back- 
ground reactions to the signal have been considered [78]. This analysis showed that the 
inclusion of background processes reduces the number of signal events by a factor of 10. 
The authors in [78] find that an e+e" linear coUider with c.m.s. energy a/s — 0.5 TeV will 
be sensitive to values of s'^ — \a\v/mN ~ 0.7 x 10"^. This amounts to the same level of 
sensitivity to the parameter \a\, for RrL scenarios with mjv = 250 GeV. The sensitivity 
to s^j^ could be improved by a factor of 3, i.e. s"]^ ~ 0.2 x 10~^, in proposed upgraded e'^e~ 
accelerators such as CLIC. 

A similar analysis should be envisaged to hold for future fi'^fi~ colliders, leading to 
similar findings for s^** = \b\v/mi\i- In general, we expect that the ratio of the two produc- 
tion cross sections of N2^3 at the two colliders under identical conditions of c.m.s. energy 
and luminosity will give a direct measure of the ratio of |ap/|6p. This information, to- 
gether with that obtained from low-energy LFV observablcs, Ou f3 f3-dccaj experiments, and 
neutrino data, will significantly constrain the parameters of the RrL models. Finally, since 
the heavy Majorana neutrinos A^2,3 play an important syncrgctic role in resonantly enhanc- 
ing potentially large CP asymmetries in their decays will determine the theoretical 
parameters of these models further. Evidently, more detailed studies are needed before one 
could reach a definite conclusion concerning the exciting possibility that electrowcak-scale 
RrL models may naturally constitute a laboratory testable solution to the cosmological 
problem of the BAU. 

7 Conclusions 

We have studied a novel variant of RL, which may take place at the electroweak phase 
transition. This RL variant gives rise to a number of phenomenologically testable signatures 
for low-energy experiments and future high-energy colliders. The new RL scenario under 
study makes use of the property that, in addition to B — L number, sphalerons preserve 
the individual quantum numbers |S — I/e,/i,r [28]. The observed BAU can be produced by 
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lepton-to-baryon conversion of an individual lepton number. For the case of the r-lepton 
number this mechanism has been called resonant r-leptogenesis [26]. 

In studying leptogenesis, we have extended previous analyses of the relevant network 
of BEs. More explicitly, we have consistently taken into account SM chemical potential 
effects, as well as effects from out of equilibrium sphalerons and single lepton flavours. 
In particular, we have found that single lepton flavour effects become very important in 
RtL models. In this case, the difference between our improved formalism of BEs and 
the usual formalism followed in the literature could be dramatic. The predictions of the 
usual formalism could lead to an erroneous result which is suppressed by many orders of 
magnitude. The suppression factor could be enormous of order 10^^ for the RrL scenarios 
considered in Section 5. Even within leptogenesis models with a mild hierarchy between the 
heavy neutrino masses, the usual formalism turns out to be inadequate to properly treat 
single lepton flavour efl^ccts; its predictions may differ even up to one order of magnitude 
with respect to those obtained with our improved formalism. 

One generic feature of RrL models is that their predictions for the final baryon 
asymmetry are almost independent of the initial values for the primordial S-number, L- 
number and heavy Majorana neutrino abundances. Specifically, we have investigated the 
dependence of the BAU on the initial conditions, as a function of the heavy neutrino mass 
scale rriN. We have found that for mjv ~ 250 GeV, the dependence of the BAU is always less 
than 7%, even if the initial baryon asymmetry is as large as 77^ = 10~^ at z — ruN/T — 0.1. 
For smaller values of itin, this dependence starts getting larger. Thus, for = 100 GeV, 
the dependence of the final baryon asymmetry on the initial conditions is stronger, unless 
the primordial baryon asymmetry is smaller than ~ 10~^ at z — 0.1. 

In order to have successful leptogenesis in the RrL models under study, the heavy 
Majorana neutrinos are required to be nearly degenerate. This nearly degenerate heavy 
neutrino mass spectrum may be obtained by enforcing an SO (3) symmetry, which is ex- 
plicitly broken by the Yukawa interactions to a particular SO (2) sub-group isomorphic to 
a lepton- type group 11(1);. The approximate breaking of U(l)i, which could result from 
a FN mechanism, leads to a Yukawa texture that accounts for the existing neutrino os- 
cillation data, except those from the LSND experiment [79]. Our choice of the breaking 
parameters was motivated by the naturalness of the light and heavy neutrino sectors. To 
obtain natural RrL models, we have followed the principle that there should be no exces- 
sive cancellations between tree-level and radiative or thermal effects. In this way, we have 
found that RrL models strongly favour a light neutrino mass spectrum with an inverted 
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hierarchy. Moreover, when the same naturalness condition is apphed to the heavy neutrino 
sector, a particular hierarchy for the mass differences of the heavy Majorana neutrinos is 
obtained. In particular, the mass difference of one pair of heavy Majorana neutrinos is 
much smaller than the other two possible pairs. 

RrL models offer a number of testable phenomenological signatures for low-energy ex- 
periments and future high-energy colliders. These models contain electroweak- scale heavy 
Majorana neutrinos with appreciable couplings to electrons and muons, e.g. Ni^2- Specif- 
ically, the (normalized to the SM) H^^-boson couplings of electrons and muons to the 
heavy Majorana neutrinos A^i,2 could be as large as 0.01, for mjv^ ^ = 100-300 GeV. As a 
consequence, these heavy Majorana particles can be produced at future e~^e~ and 
colliders, operating with a c.m.s. energy ^/s = 0.5-1 TeV. Another feature of RrL models 
is that thanks to the inverted hierarchic structure of the light neutrino mass spectrum, 
they can account for sizeable Ou^f] decay. The predicted effective neutrino mass |(m'')ee| 
can be as large as 0.02 eV, which is within the sensitivity of the proposed next round 
of 0i//9/9 decay experiments. The most striking phenomenological feature of 3-generation 
(non-supersymmetric) RrL models is that they can predict e- and /x-number-violating pro- 
cesses, such as the decay fi ^ e'j and fi —>■ e conversion in nuclei, with observable rates. 
In particular, these LFV effects could be as large as 10~^^ for B{ii ej) and as large as 
5 X 10~^^ for a — > e conversion rate in 22Ti, normalized to the /i capture rate. The above 
predicted values are within reach of the experiments proposed by the MEG and MECO 
collaborations. 

Although the present study improves previous analyses of the BEs related to leptoge- 
nesis models, there are still some additional smaller but relevant effects that would require 
special treatment. The first obvious improvement would be to calculate the thermal effects 
on the collision terms, beyond the HTL approximation. These corrections would eliminate 
some of the uncertainties pertinent to the actual choice of the IR regulator in some of 
the collision terms. These effects limit the accuracy of our predictions and introduce an 
estimated theoretical uncertainty of 30% for leptogenesis models operating well above the 
electroweak phase transition, with relatively large K factors, i.e. Ki^^ > 5. For models at 
the electroweak phase transition, the IR problem is less serious, but larger uncertainties 
may enter due to the lack of a satisfactorily accurate quantitative framework for sphaleron 
dynamics. Although the implementation of the sphaleron dynamics in our BEs for RL 
models was based on the calculations of [27,30,33], particular treatment would be needed, 
if the electroweak phase transition was a strong first-order one. In this case, the dynamics 
of the expanding bubbles during the electroweak phase transition becomes relevant [80]. 
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This possibility may emerge in supersymmetric versions of RL models. Nevertheless, the 
inclusion of the aforementioned additional effects is expected not to modify the main results 
of the present analysis drastically and will be the subject of a future communication. 
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A Collision Terms 



A.l Useful Notation and Definitions 

The following notation and definitions are used in the derivation of the BEs. The number 
density, Ua, of a particle species, a, with ga internal degrees of freedom is given by [45] 

n„(T) = I ^ exp [ - {^Jp' + ml - l^a{T))/T 

gamlTe^'-^^y^ /ma\ ^^^^ 



27r2 " V T 

where /la is the T-dependent chemical potential and Kn{x) is the nth-order modified Bessel 
function [81]. In our minimal leptogenesis model, the factors ga are: gw^ — ^93 = 6 and 
9^ = 5'$t = 2, and for the ith family: gN^ = 2, gi, = g^c = 4, gg. = ggc = 12, and 
9ui — 9u'? — 6. Using the same formalism as [17] the CP-conserving collision term for a 
generic process X ^ Y and its CP-conjugate X ^ Y is defined as 

7# = 7(X ^Y) + 7(X ^ F) , (A.2) 

with 

7(X ^ y) = ldnxd7TY{27T)U^^\px-PY) e~Px/T \M{X ^Y)\\ (A.3) 

In the above, |7V1(X y)p is the squared matrix element which is summed but not 
averaged over the internal degrees of freedom of the initial and final multiparticle states X 
and Y. Moreover, dnx represents the phase space factor of a multiparticle state X, 

where Sx — JT-id' is a symmetry factor depending on the number of identical particles, nid, 
contained in X. 

As CPT is preserved, the CP-conserving collision term 7y obeys the relation 

lY -ll- (A.5) 
Analogously, it is possible to define a CP-violating collision term ^7^ as 

57# ^ liX^Y) - 7(X^F) = -57]^, (A.6) 
where the last equality follows from CPT invariance. 
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A. 2 CP-Conserving Collision Terms 

In numerically solving the BEs, we introduce the dimensionless parameters: 



^ = -7^ , X = , = , Qr ^ [ , C„ = , A.7) 



where a — 1,2,3 labels the heavy Majorana neutrino states, s is the usual Mandelstam 
variable and miR, is an infra-red (IR) mass regulator which is discussed below. 

In terms of the resummed effective Yukawa couplings (h'^)ia introduced in [17], the 
radiatively corrected decay width of a heavy Majorana neutrino into a lepton 
flavour I is given by 

= ^ liK)iiK)ia + mimia] . (A.8) 

By means of (A. 3), the 1 —>■ 2 CP-conserving collision term 7^1, is found to be 

4 / J 

= ^^^'^^V^- K,{zV^), (A.9) 



TT^ Z 



where E^^^p) = yp^ + m^^ and g^^ = 2 is the number of internal degrees of freedom 
of Na- Upon summation over lepton flavours /, this collision term reduces to the corre- 
sponding one given in (B.4) of [17]. 

For 2 — > 2 processes, one can make use of the reduced cross section a{s) defined as 
a{s) = 87r$(s) j dny {2tiY 5^^\q - Py) \M{X^Y)\^ , (A.IO) 
where s — and the initial phase space integral is given by 

= Jd7rxi27r)U^'^{px-q). (A.ll) 

These expressions simplify to give 

1 *^ 

= 8^ /^^ \M{X^Y)f , (A.12) 
t- 

where t is the usual Mandelstam variable, and the phase-space integration limits t± will be 
specified below. 
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In processes, such as NaV^ Li^, the exchanged particles (e.g. L and $) occurring 
in the t and u channels are massless. These collision terms possess IR divergences at the 
phase-space integration limits t± in (A. 12). Within a more appropriate framework, such 
as finite temperature field theory, these IR singularities would have been regulated by the 
thermal masses of the exchanged particles. In our T = field theory calculation, we have 
regulated the IR divergences by cutting off the phase-space integration limits t± using a 
universal thermal regulator mjR related to the expected thermal masses of the exchanged 
particles. This procedure preserves chirality and gauge invariance, as would be expected 
within the framework of a finite temperature field theory [43]. 

Thermal masses for the Higgs and leptons are predominantly generated by gauge and 
top-quark Yukawa interactions. In the HTL approximation, they are given by [44] 

mjiT) 1/2 /2\ 



2^2 

ml{T) 
2^2 




(A.13) 



where d — (8M^ -|- M| -|- 2mf -|- Mfj)/{%v'^). In our numerical estimates, we choose the 
regulator mjR to vary between the lepton and Higgs thermal masses, evaluated at T m^v- 
The resulting variation in the predicted baryon asymmetry can be taken as a contribution 
to the theoretical uncertainties in our zero temperature calculation. 

For reduced cross-sections with an apparent singularity at the upper limit the 
following upper and lower limits are used: 

t+ = -m^R, = m%^-s. (A. 14) 

Likewise, for reduced cross-sections with apparent singularities at both the upper and lower 
limits t±, the following limits are employed: 

It is important to remark here that the collision terms do not suffer from IR singularities 
at T < Tc, because the leptons, W and Z bosons receive f(T) -dependent masses during 
the electroweak phase transition. The full implementation of such effects will be given 
elsewhere. 

Substituting (A. 10) and (A. 11) into (A. 3), one obtains 

4 oo 

= / ^ ^i(^V^) ^rix) , (A.16) 

d4 n z J 
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where Xthr is the kinematic threshold for a given 2 — > 2 process. 

For 2^2 AL — 1 processes, one can repeat the procedure in [17] (Appendix B), 
without summing over lepton flavours. Each AL = 1 process has an identical factor 
dependent on h'^. To produce the AL — 1 collision terms for each lepton flavour, this 
factor needs to be replaced with its un-summed equivalent, 

{K)l(K)l'^ + {h^-TlaOaia. (A.17) 

exactly as was done in (A. 8). The remainder of the analytic expression for each of these 
terms is presented in [17]. 

In addition to the Higgs and gauge mediated AL = 1 terms, there arc also 2^2 
AL = 2 processes. As before, these processes are Lk^ ^ -^F"^^ and LkLi ^ <|)t$t where 
the former has its real intermediate states subtracted. The analytic forms of these collision 
terms are identical to the total lepton number case but lepton flavour is not summed over. 
The reduced cross sections are given by 



E Re 

a,/3=l 



+ 2 



k/3 



(tt) 



and 



$t#t 



E Re (hll hip h\* h'{pj Bafi , 
where the A and B factors are presented in [17]. 



(A.18) 



(A.19) 



As we now consider lepton flavours separately, it is necessary to include AL = 0, but 
lepton flavour violating interactions. The three lowest order 2 <-> 2 processes are shown 
diagrammatically in Figure 10: Lfc$ ^ Li^, Lfe$^ ^ L;$^ and L^Lp ^ (note that 
k ^ I). The flrst of these reactions contains heavy Majorana neutrinos as RISs. These 
need be removed using the procedure outlined in [17]. The reduced cross section for each 
of these processes is 



Li* 



E [ih 

a,/3=l 



la 



Ch-)la iKh Ch-)kp + Ch-Tla iK)la (h-h iKU ^a/J (A.20) 



with 



Xttr 



Call 



(« = /3) 



(A.21) 
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--$t Lfe 



Li L\ 



c 



* $ 



(a) 



(b) 



(c) 



Figure 10: AL = interactions between leptons of different flavours mediated by heavy 
Majorana neutrinos, (a) corresponds to the process L^^ <-> L;$, (b) corresponds to the 
process L^^^ <-> L;$^ and (c) corresponds to the process L^Lf <-> 



In (A. 21), Pq, ^{x) is the Breit-Wigner s-channel propagator 

1 



(A.22) 



Therefore, foUowing the procedure in [17], the RIS-subtracted propagator is determined by 



TT 



(A.23) 



Processes (b) and (c) in Fig. 10 do not contain RISs and have the following reduced 
cross sections: 

3 



a 



$t$ 



3 



a,/3=l 



a/3 ) 



(A.24) 
(A.25) 



where ior a ^ (5, 



V. 



SnR — 



■Kx{aa - ap) 



( x + aA (x + ac 

[x + a/3) In I — ] - [X + Qa) In 



In 



and for a — (3, 



X 



f — 

7r(x + aa) 



(A.26) 
(A.27) 

(A.28) 
(A.29) 
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